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Abstract 

Motivated by an analogy with the spin anisotropies in the quantum XY chain and 
its reformulation in terms of spin-less Majorana fermions, its bosonic analogue, the spin- 
anisotropic quantum spherical model, is introduced. The exact solution of the model 
permits to analyse the influence of the spin-anisotropy on the phase diagram and the 
universality of the critical behaviour in a new way, since the interactions of the quantum 
spins and their conjugate momenta create new effects. At zero temperature, a quantum 
critical line is found, which is in the same universality class as the thermal phase transition 
in the classical spherical model in d -|- 1 dimensions. The location of this quantum critical 
line shows a re-entrant quantum phase transition for dimensions 1 < d < 2.065. 
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1 Introduction 


The study of equilibrium phase transitions has taken enormous benefits from the analysis of 
exactly solvable models [l5l [TJ |39l [62l [TTl [13]. The classical spherical model, invented in 
the seminal work by Berlin and Kac [8] and with its subsequent simplification by Lewis and 
Wannier [52], has been a valuable test system for the explicit analytical verification of more 
general scaling descriptions, in a specific setting (examples are critical behaviour of observables 
or finite-size scaling). It is related with more realistic spin systems as the n —)■ cxo limit of the 
0(n)-symmetric Heisenberg model [66]. It is well-known, as already observed by Berlin and 
Kac [8], that in the original formulation in terms of classical spin variables Si G M, the specific 
heat does not vanish in the zero-temperature limit, hence the Nernst theorem is not obeyed 
in this model. This was one motivation to take the quantum nature of the spin variables into 
account, by a canonical quantisation scheme, and has lead to Obermair’s formulation of the 
quantum spherical model m- This takes the form of a quantum rotor model, where the kinetic 
energy term in the hamiltonian does not commute with the spin-exchange interactions. The 
properties of this exactly solvable model have been analysed in great detail, see e.g. [SHElllMl 
EHiiiHiEa nniaaisHiini [To]. Independently, the quantum spherical model was also obtained 
via the so-called ‘hamiltonian limit’ [50] as the logarithm of the transfer matrix in an extremely 
anisotropic limit [6^ [3^ [38] . In its most conventional formulation as a quantum rotor model 
[571 [35l [68] , the quantum spherical model may be obtained as the limit n —)■ cxo of the quantum 
non-linear 0(n) sigma-model [68]. For different choices of the kinetic energy term, there are 
further quantum spherical models, which become the n —)■ cxo limit of an SU(?7,) Heisenberg 
ferromagnet or anti-ferromagnet [5B1132] . 

Quantum spherical models have been discussed in the context of specific applications, for 
example for the description of networks of Josephson-junction arrays [27] E]. Certain modern 
theories of cuprate supraconductivity are based on SO(5)-symmetric quantum non-linear sigma 
models, and it is thought that this kind of models might be an effective description of the large- 
distance, low-energy properties of more realistic models, see e.g. [22] for a detailed review. 

Habitually, (mean) spherical models are defined in terms of a classical hamiltonian 


T^ci = 
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-jYiSr 

i=i 




n^n-\-e 


- BSr, + 


( 1 . 1 ) 


with the spherical spins Sn G M. Herein, n runs over the sites of a d-dimensional hyper-cubic 
lattice with J\f = N'^ sites, the vectors e^, 1 < j < d, are the unit vectors in the direction, 
B is an external magnetic field and J is the exchange integral. Finally, the spherical spins obey 
the mean ‘spherical constraint’ {S^) = AT [H] [S2], from which fi is found. 

The generalisation towards a quantum spherical model is formulated by considering now the 
spins Sn Sn as operators, and introducing canonically conjugate momenta Pn, which obey 
the canonical commutation relations 
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( 1 . 2 ) 


The most common ansatz for the quantum hamiltonian is to make "Hci ^ Hd an operator and 
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to add a kinetic energy term of non-interacting momentaji] viz. 






n+e, 


BS„ + + |p2 


(1,3) 


with a new coupling g, which controls the strength of the quantum fluctuations. In equilibrium, 
one can express the spherical constrainll^ as a thermodynamic derivative 



2 dlnZ 
T dn 



(1.4) 


where Z = trexp(—iJ/T) is the partition function and T is the temperature. This quan¬ 
tum hamiltonian can also be obtained as the logarithm of the transfer matrix of the classical 
spherical model in d -|- 1 dimensions, in a certain strongly anisotropic limit [23 EHl ESI E21 E2]- 
This mapping in particular shows that the zero-temperature quantum critical behaviour of the 
quantum phase transition of the ground-state of the quantum spherical model fll.Sp in d dimen¬ 
sions [35] is in the same universality class as the hnite-temperature transition of the classical 
spherical model in d -|- 1 dimensions [53 ESI E21 iHl ESI ESj- 

It is straightforward to recast the hamiltonian fll.3p in terms of bosonic ladder operators an 
and ajj, dehned as follows 



which obey the canonical commutator relations 

^n.m i \^ni ®m] \p\ii ^ 

and render the hamiltonian fll.31) as follows 

r (alan + 0 - {an + a^) 

d 

J=1 



(1.5) 


( 1 . 6 ) 


(1.7) 


The computation of the eigenvalues of such hamiltonians is a matter of finding the appropriate 
canonical transformation and is treated in appendix A. Here, we wish to point out an analogy 
with quantum Ising/XY chains (also called Ising/XY chains in a transverse field), with an 

^Even in the case of competing interactions, where new multicritical points, called Lifshitz points [43], 
can be found in the classical spherical model and which may present strongly anisotropic scaling behaviour, see 
[29l|30l|23liMlmuss] and refs, therein, existing studies on the quantum version do not consider any interactions 
between the momenta m- 

^Sometimes the constraint is given in the form J2n “ A”/!, see e.g. [68l [58], which in the zero 

temperature limit amounts essentially to a re-scaling of the spherical parameter. Throughout, units are such 
that the Boltzmann constant fcs = 1. 
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anisotropy in spin space, and given by the hamiltonian | 1 S 1 E] 
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( 1 , 8 ) 

(1.9) 


where the denote the Panli matrices attached to the site of a periodic chain of N sites. 
The transverse field g measures the quantum fluctuations and A is a spin-anisotropy coupling. 
After a Jordan-Wigner transformation, the hamiltonian (11.81) can be brought to a quadratic 
form fll.91) in the fermionic ladder operators c„ and (we did not carefully specify the non-local 
boundary conditions in the fermionic variables since we shall not require their form) with the 
anticommutator relations 


{Cn,Hn} = ^n,m , {c„, C^} = ^} = 0 (1.10) 

The ground-state of quantum Ising/XY chain fll.Sp has a rich phase diagram with a disordered 
phase for > 1 , a line of second-order transitions at 5 ^ = 1 which is in the universality class 
of the 2D Ising model for A 7 ^ 0, an ordered ferromagnetic phase for \/l — A^ < < 1 and an 

ordered oscillating phase for g < \/l — A^ [ 6 l (STJ [T^ EHl El 125] . The universality of the quantum 
critical behaviour at T = 0 , including the universal amplitude combinations [59l [ 6 OI 001 ITT] , 
with respect to 0 < A < 1 along the Ising critical line has been explicitly conhrmed: for the 
chain for both the spin-A as well as the the spin-1 representations of the Lie algebra of the 
rotation group [37], as well as in 2D for the spin-A representation [36] . 

Comparing the fermionic hamiltonian (II.Oh with the bosonic one ( 01)0 one observes that 
in the former the two-particle annihilation/creation processes are controlled by the parameter 
A, whereas that parameter happens to be fixed to unity in the latter. Here, we shall inquire 
into what happens if an analogous rate is introduced into the hamiltonian fll.7p . and write 


( 1 . 11 ) 

A)fi'Pn^n,+ej) 

The re-formulation in terms of the original spins and momenta shows that the hamiltonian 
fll.lip introduces an interaction between the momenta, quite analogous to the spin anisotropies 
in the quantum XY chain fll.Sp . In the special case A = 1, this new interaction disappears and 
one is back to the quantum rotor spherical model as studied in the literature so far. We call 
the model defined by (11.111) the spin-anisotropic quantum spherical model (saqsm), because of 
the analogy of the parameter A with the spin anisotropy in the fermionic hamiltonian (II.8111.9p . 

^Alternatively, one can consider the fermionic degrees of freedom in (HH) as hard-core bosons. Relaxing the 

‘hard-core/fermionic’ constraint on the single-site occupation numbers (hi) = (cjci) = 0,1, towards = 

vAf, where i/ = i is a filling factor, one has a third way to replace (HH) by a quantum spherical model |55] . 
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( 1 . 12 ) 


It will be convenient to work with the spherical parameter (already used in fll.lip i 



For B = 0, there is a duality transformation Sn t-)- Pn, ^ g, X ^ —A. It is therefore sufficient 
to restrict attention to the case A > 0, as we shall do from now on. In the special case A = 0, 
pairs of particles can neither be created, nor destroyed, which formally is expressed through 
the conservation, expressed by [N,H] = 0, of the total number of particles N := 

This case has properties different from the situation where A 7 ^ 00 

This work is organised as follows. Section 2 presents the general formalism for the solution 
of the model and the new techniques required for its analysis when A 7 ^ 0,1. We shall focus 
on the quantum phase transition at zero temperature. A detailed analysis of the spherical 
constraint surprisingly shows that for dimensions 1 < d < 2.065, there is a re-entrant quantum 
phase transitions when A is small enough. There is no known classical analogue of this effect. 
The critical behaviour and its universality along the A-dependent critical lines will be analysed 
and we shall discuss the relationship with the thermal phase transition of the classical spherical 
model. As one should have expected, we hnd a critical line 0 for 0 < A < 1, where the quantum 
critical behaviour of the SAQSM is in the same universality class as in the classical spherical 
model in d -I- 1 dimensions. Section 3 gives our conclusions. Technical details are treated in 
several appendices. Appendix A recalls the exact diagonalisation techniques, in appendices B 
and C the spherical constraint and the consequences for the quantum critical point are studied, 
in appendix D the spin-spin correlator is derived and appendix E looks in more detail into the 
existence of the re-entrant quantum phase transition. 


2 Solution and quantum phase transition 

2.1 General formalism 


In order to analyse the thermodynamic behaviour of the quantum spherical model fll.lip with 
A arbitrary, the first task is to bring H into a diagonal form. This calculation is carried out in 
appendix A, and leads to 

H = ^2h^gJ/sJ2^k (bih + l)+Ho (2.1) 

k&K ^ ' 

where the eigenvalues are given in eq. (IA.9p 


At, . sl\.u 


\ 


1 + X 

s - — 2_^ cos kj 

i=i 


\ 


1 — A 

s- — 2_^ cos kj 

i=i 


and the quasi-momenta 1C 3 kj = with rij = 0,1,... A^ — 1 and j = 1,. 
reciprocal lattice /C. Finally, from flA.lGp we have 


Hn = 


U s - {1 + X)d/{2) 


( 2 . 2 ) 

., d, with the 
(2.3) 


"‘The conservation of N is reminiscent of the spherical constraints used in [56l|32], although the quantum 
critical behaviour of the A = 0 model (11.111) will turn out to be different. 

®Our methods of analysis are restricted to |A| < 1, see appendices B and C. 
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Since the quasi-particles are independent, non-interacting particles, the calculation of the par¬ 
tition function reduces to a computation of products of geometric series, such that the free 
energy F = —T\nZ reads explicitly 


F 


TM\n2 - 


4Js- (1 +A)d/2 


-|- T In sinh 

h gJ 
TV 2s \ 

1 -I- A 

s o Z^cos/cj 

1 — A 

s o 

k 


i=i \ 

3 = 1 


(2.4) 


At this point, one can go to the inhnite-size limit M = N'^ ^ oo. In particular, the spherical 
constraint fll.dp then takes the form 



dk 

{2tiY 


coth 



s — 


i-y 

As 


Ell cos kj 


2Ai 



(2.5) 


where B = [—tt, tt]'^ is the Brillouin zone. For the special case A = 1, we recover the form of the 
spherical constraint known from the literature, see [351I6HII6211IIIEH]. 


Besides thermodynamic observables, we shall also study the spin-spin correlator. In ap¬ 
pendix D, it is shown that 


{SnSn+r) — 


n?g 


dk 


2s — (1 — A) Eli cos k 


8Js Js {2nY ^ 2s - (1 + A) E,li cos k, 


coth 


l2hVVsAfc/(2T)J JJcos(rj/cj) 

i=i 

( 2 . 6 ) 


2.2 Quantum phase transition 

In d > 2 dimensions, the spherical model undergoes a phase transition at some critical tem¬ 
perature Tc > 0 [68l [55l EH [m EH EH]- In general, one expects that this hnite-temperature 
transition of the d-dimensional model should be in the same universality class as the one of 
the classical model (without quantum terms) [SDIEHIIII- Here, we rather concentrate on the 
quantum phase transition which occurs in the ground-state, that is, at temperature T = 0. 

Generically, quantum phase transitions arise mathematically from a degeneracy in the 
ground-state of the hamiltonian. In order to localise the quantum critical point in terms of the 
model’s parameters, consider the smallest energy gap AF' 


AF' := lim A^ = 

k^O 



(2.7) 


This energy gap closes for 



( 2 . 8 ) 


such that the spherical parameter must satisfy s > (1-|- |A|)d/2. The ground-state thermody¬ 
namics now follows from an analysis of the spherical constraint (I2.5p . which in the limit T —)■ 0 
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takes the form 



dk 


2 _ 1-A^ 

^ A 


Ell cos kj 


+ 


¥ Ell cos fc,] V2 [, - J2U 11^^ 


£ 1 "l 

2Js-i±Arfy 
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1 (2.9) 


This dehnes the function s = s{g, X,d, B), or alternatively its inverse g = g{s, X,d, B). For a 
vanishing external held 5 = 0, this equation is symmetric under A —A, hence it is then 
sufficient to consider the case A > 0 only. We shall almost always restrict to this special case, 
and then write g = g{s, A, d) := g{s, A, d, 0). 

1. For A = 0, the constraint simplihes considerably and can be worked out explicitly 


1 = 


B 


2Js-d/2 


gh? 


dk 


8Js^ Jq {27ry 


+ ^J^cosfcj 


i=i 


B 


2Js-d/2 


+ 


8Js 


( 2 . 10 ) 


Eq. fl2.10p gives directly the inverse function g = g{s, 0, d), where d appears as a real parameter. 

2. For A = 1, this has been analysed many times and it is well-known [3^ l 68 t ITSt [TH [58] 
that fl2.9p can be re-written as (set 5 = 0) 



2\d 


du /o(W) 


( 2 , 11 ) 


where Iq is a modihed Bessel function [T]. Again, this formulation has the appealing feature that 
by now d can be considered as a continuous parameter in an analytic continuation g = g{s,l,d). 

3. Finally, for generic A, the constraint ( 12 . 9 p can be written in the form (set 5 = 0) 


/SttV 

n^g 


= S 2 


du 




'0 


X 


— d{d — 1 ) 


\/x{l - x) 

l-X^higf dl-X^ 


1 + 


h{0) 


4 Jo(p)2 2 4 V Io{g) 

where the /„ are modihed Bessel functions [1] and we dehned the function 

/ \\ /l-|-A , ,1 — A 

g = g{u, x. A) := M I x — -H (1 — x )—-— 


( 2 . 12 ) 


(2.13) 


Eqs. (12.12112.13p contain d as a real parameter and give directly g = g{s, A, d). This form of the 
constraint is derived in appendix B. 


2.3 Critical behaviour 


Now, the constraints fl2.10p . (12.111) and fl2.12p can be used to extract the quantum critical 
coupling and the relation between g and the spherical parameter s for the diherent values of A. 

1. First, we consider the case A = 0. From fl2.10p . we have, even for 5 7 ^: 0 



(2,14) 
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With the critical value = d/2, see eq. fl2.8p . we have the critical coupling, for i? = 0 


gc = gc{^,d) ■.= g{sc,Q,d) = (2.15) 

which is non-vanishing for any dimension d > 0. For the later extraction of the critical expo¬ 
nents, we also note {g — gc)/gc = (2/d)(s —Sc)- This linear behaviour is independent of d, hence 
there is no upper critical dimension. 

2. Next, we briefly recall the known result for A = 1. We are interested in hnding the 
critical value gc = 5'c(l, d) := g{sc, 1, d), if it exists and to obtain the variation of g close to gc, 
which we can describe in terms of 


, _ (— _ [^9c^ 

" ■ V V ^3/2 


(2.16) 


Consider a := s — Sc = s — d. In order to extract from fl2.11l) any non-analytic terms in a, 
one may formally split [55] the domain of integration The first term, if it 

exists, will give a analytic contribution to g{s) near s ~ Sc, in particular g^, = g{sc) in the limit 
r] —)■ oo] the second term will give any non-analytic contributions which may arise. In order to 
hnd those, recall the asymptotic form /o(p) — e^(27rp)“^/^ as p ^ cxo [1]. Then, for d < 3 


^ (2v/F)-^ ^ r due-- 1. 

® (27r)'^/2 (27r)C+i)/2y2 yC-i)/2+i ('27]-)C+i)/2 \ 2 J ^/2 

(2.17) 

(the Gamma function r(x) [1] is dehned via analytic continuation, if needed) and only now 
one also lets rj —)■ oo. For d > 3, tg ^ a + is dominated by the analytic term. 

Finally, for d = 3, the non-integrability gives rise to a logarithmic correction such that hnally 

[MlEHlIIHlEainiEHliniilQ] 


( ; ifd<3 

-—— ~ fg ~ < ^3 (T In cr ; if d = 3 
y Ay a ; if d > 3 

as cr —)■ 0 and with known constant amplitudes A^, A^, Ay, see appendix C. 
Explicitly, the critical coupling gc{l,d) can be expressed as an integral 


(2.18) 


Pc(l,d) = 271 



J 

¥ 


(2.19) 


The asymptotic behaviour of /o(p) for P large tells us that pc(l, d) > 0 is hnite for d > 1 but that 
Pc(l,d) = 0 for d < 1. For d = 2, the identities [HU eq. (2.15.20.5)], P eqs. (8.1.2),(15.1.26)] 
give the closed expression 


dc(l,2) 



J J 

— ~ 9.67826 — 


( 2 . 20 ) 


which agrees with the numerical values quoted in [TS] [55] • The result fl2.20p is the counterpart 
to the exact value of Tc in the 3D classical spherical model [T5] . 
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Figure 1: Critical coupling gc{X,2) from 
ation/annihilation rate A, for d = 2 space dimensions. 


as a function of the pair cre- 


3. In the general case 0 < A < 1, the asymptotic analysis of the spherical constraint is 
more involved than in the two previous cases. As far as the critical exponents are concerned, 
we show in appendix C that 02.181) remains valid for 0 < A < 1, where the amplitudes are given 
explicitly by eqs. 0(1.130(1.140(1.151) FI 


Turning to the values of the critical coupling = gc{X,d), we consider hrst the 2D case 
and have 

r?(A,2) = 87r2&^^4 (2-21) 


where, using 
G(A) = / 


f du [ 

0 Jo 

rl 


G(A)2 

eq. (2.15.20.5)], we hnd from 02.12p 
dxexp (—uSc) 




dx 


lo 

1 - A 
16 

1-A 


X) 


X) 

1 + 3A 


—-—(1 + 3A)/o(p)^--— Ii{qY H--A(^?)-^o(^?) 


2p 


1 — X 


2A 

ITa 




2A 


3-^2 


,l,-;2,2: l-a;- 
2’ ’2’ ’ ’ V 1 + A. 


( 2 . 22 ) 


This quite explicit form is more easily treated numerically than the full double integral 02.12p . 


®For A > 1, the asymptotic methods used in appendix C for analysing (12.121) cannot be taken over, since the 
argument g, see eq. (12.131) . of the Bessel functions can vanish. The contributions of such zeroes would have to 
be included into the analysis. However, since the numerical values do not show evidence for a singularity at 
A = 1, we expect that our results should be straightforwardly generalisable to A > 1. 





































gc{>‘,d) [Jh 2 




Figure 2: Left panel: Critical coupling gdX, d), computed from fl2.23l) . as a function of A for d = 
[1.3,1.4,1.5,1.6,1.7,1.8,1.9, 2.0, 2.1] from bottom to top. Right panel: slope dg^X, d)/d\\^^Q 
of the critical coupling at A = 0, as a function of d. For d ~ 2.065, the slope vanishes. 


to be considered in generic dimensions d. In hgure [U, we plot gc = gdX, 2) over against A. 
While the two known values fl2.15112.201) for A = 0 and A = 1 are certainly reproduced, we 
also observe that the behaviour of gdX, 2) is not monotonous in A, but rather has a minimum 
around A ~ 0.1. This surprising feature of a re-entrant quantum phase transition does not have 
an analogue in the 3D classical spherical model. 

Indeed, this re-entrant transition for A small enough is a generic feature of the quantum 
spherical model. In the left panel of hgure [2], we show the critical coupling gdX, d), as given by 


5'c(A,d) 



exp(—n(l -|- X)d/2) 
\/x{l - x) 


X 


d{d — 1) 


4 loiQd 


dl-X^ 
2 4 



Io{qY 

h{g) \ 



(2.23) 


Clearly, the hgure suggests that gdX, d) should go through a non-vanishing minimum for all 
dimensions d ;< 2.1. 

Let us make this statement more precise. First, we observe from fl2.20p that gd^,2.) > 
gd^, 2) = 8J/h?. Second, from fl2.19p it follows that gd^, d) grows monotonously with d. Since 
gdX, 2) is increasing with A for A large enough, see hgure [H this means that the slope of gdX, d) 
at A = 1 should be positive, viz. dgdX, d)/dX\y^^^ > 0. On the other hand, in appendix E we 
show that close to A = 0 one has 


9d\d) ~ gd^^d) 


5 f(o)A'^/^ ; if 1 < d < 2 

5 f(i)A ; if d > 2 


(2.24) 


and where the known constant ^'(o) > 0, but the sign of the known constant g(i) may depend on 
d. Therefore, the slope dgdX,d)/dX\y^^Q < 0 for dimensions 1 < d < 2 and diverges as A —)■ 0. 
On the other hand, in the right panel of hgure [21 we show the hnite slope dgdX, d)/dX\^^Q of gc 
at A = 0, for dimensions d > 2, as a function of d. Clearly, the slope oi gc at X = 0 is negative 
for d ;< 2.065 and becomes positive for larger values of d. For d small enough, the slope of 
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gc{X,d) is negative at A = 0 and positive at A = 1. By Rolle’s theorem, the critical coupling 
5 'c(A, d) should have a minimum at some non-vanishing value of A, for all dimensions d < 2.065. 
This is indeed what we observe in the left panel of figure [2l In consequence, the spin-anisotropic 
quantum spherical model has a re-entrant quantum phase transition for dimensions d < 2.065. 


2.4 Physical observables near quantum criticality 


The scaling of the thermodynamic observables follows from the free-energy density. Since we 
restrict ourselves to an analysis of the zero-temperature properties of our model, the quantum 
coupling g takes over the role of the temperature in classical spin systems, such that tg as 
dehned in fl2.16112.ISp takes over the role of T — in classical phase transitions. Therefore, one 
expects for the singular part of the free energy density 


/ = 


F 

Jf 


+ W 2l 




1 + A 
2 


d 

cos kj 
i=i 



1 - A 
2 


E 

i=i 


cos kj 


(2.25) 


to obey the following scaling behaviour 


nitg,B) = A^\tg\^-'^Wi {A2B\tg\-f^-^) 


(2.26) 


where W± are universal scaling functions, associated with the sign of tg ^ 0 ) and a, /3 ,7 are 
the standard critical exponents. All non-universal information on the specihc model can be 
absorbed into the two metric factors Ai 2 . Similarly, we consider the spin-spin correlation fl2.6p 
^(l’"l) = (SnSn-er) at zero temperature T = 0. As shown in appendix D, we can use spatial 
translation- and rotation-invariance, and have for A > 0 


C{R) = (SoSr) = 


Js 


dk 2s - (1 - A) Y.'j=i cos kj 
2s-(1 + A) 


cos kiR 


h?g 1 s - ^d f 1 5^ 

7^(27r)('^+b/2^A(l + A)rf/2 


K.-. I I 


where we identify the correlation length, with s = |(1 + X)d + a, as follows 

'iTX 




a 


-1/2 


(2.27) 


(2.28) 


and Ky{x) is the other modihed Bessel function [T]. For isotropic classical phase transitions, a 
long-standing result of Privman and Fisher [59] states that there exist only two independent 
non-universal metric factors, such as Ai 2 . For quantum systems, anisotropies are possible be¬ 
tween correlators along the spatial lattice and correlations in the (euclidean) ‘time’ direction 
and generated via the transfer matrix T = exp {—tH). One then must distinguish ‘parallel’ dis¬ 
tances ry along the ‘time’ direction and ‘perpendicular’ distances r± along the space direction. 
The correlation length = .^_l considered here is spatial, whereas the ‘temporal’ correlation 
length ^11 (AE) ^ is related to the energy gap of E[. The anisotropy between ‘time’ and 
‘space’ introduces a further metric factor which in those cases where there is a classical ana¬ 
logue, and therefore the dynamical exponent z = 1, amounts simply to a further independent 
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amplitude Dq related to the freedom of normalisation of the quantum hamiltonian H. For 
such anisotropic or quantum systems (at T = 0), one expects a scaling form for a two-point 
correlator [10l[iail9] 

CiR-,tg,B) = (2.29) 

where in the situation under study here, we have R = |i?| = |rj_| and ry = 0. As before, X± 
are universal scaling functions with non-universal metric factors F*o,i, 2 - For isotropic systems, 
one has z = 1 such that the distinction between the scaling of and ry is no longer necessary 
and Dq = 1 without restriction to the generality. Then, in that situation, only two of the 
four metric factors Ai^ 2 , D 12 are independent, according to the long-standing Privman-Fisher 
hypothesis [59]. This follows by tracing the metric factors as they occur in the thermodynamic 
observables and using the static fluctuation-dissipation theorem. For potentially anisotropic 
or quantum systems, even if z = 1 , this argument has to be generalised in order to admit a 
potentially non-universal normalisation Dq. This leads to the following universal amplitude 
combinations Qi, 2,3 [ 10 ] 

Qi = Q 2 = Q 3 = (2.30) 


where the amplitude ^0 is from ^ ~ Here, we shall use the dependence on the parameter 

A > 0 to control explicitly the universality and hence to test the scaling forms fl2.26112.29p . 

Returning to the quantum spherical model at T = 0, the analysis of the spherical constraint, 
see appendix C, has given us the dependence of the shift tg on the shifted spherical parameter 
a = s — Sc- Including now the magnetic held B as well, we have to leading order in a 


t 


9 





r'^ 


cr^ ; if d < 3 

A3 a In a ; if d = 3 

a ; if d > 3 


(2.31) 


with explicitly known amplitudes A<, A 3 and A>. For a non-vanishing magnetic held B ^ 0 
the magnetic contribution will always dominate the behaviour of the spherical constraint near 
criticality. 

1. First, we treat the case 0 < A < 1 and 1 < d < 3. From the Gibbs free energy, eq. 
fl2.25p . we hnd for the magnetisation near criticality 


m{tg,B) 


df{tg,B) B 1 

dB 2J a 


(2.32) 


where the spherical constraint fl2.3ip must be used. The critical behaviour is extracted by 
moving along the quantum critical ‘isochore’ B = 0 or else the quantum critical ‘isotherm’ 
tg = 0. We obtain 


m{tg,0) ~ 




8J 


8J 


d-1 
. 5d+3 


(2.33) 


where we used the non-universal amplitudes from fl2.3ip and the value of Qc = gc{)\,d), which 
are explicitly A-dependent. 
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The analogue of the susceptibility is dehned by xi^gyB) = dm(tg, B)/dB. Explicitly, we 
hud 


x(4,o) 

x(o,5) 


2 

A d-i 2 

2J ^ 
d+3 < ^ ^ 



4 

B ^+3 


(2.34) 

(2.35) 


In general, the specihc heat is given by the second derive of the free energy with respect to 
the temperature (here replaced by tg). Here, we consider its analogue, where the role of T is 
taken over by tg. Furthermore, in the spherical model, the spherical constraint requires a little 
more careful consideration, which amounts to 


C{tg,B) 


d f dr^tg,B) \ 
dtg V dtg j 


(2.36) 


where the hrst derivative must be taken grand-canonically, with hxed spherical parameter, 
whereas the second derivative is an usual thermodynamic derivative, in the canonical ensemble, 
see e.g. [H EH m [351 HB [lo] . We End 


c(tg,0) — Co + 


c(0, B) — Co + 


d - 1 V 8J 




(1 2 d—3 


d — 1 V 8 J 


d+3 


(8J3) ^+3 


where Cq is an unimportant background constant. 

The correlation length introduced in eq. fl2.28p . reads near criticality 


(2.37) 

(2.38) 


5(ip.O) = (2-39) 

1 

5(0,-S) = (2J)ra'-B-ra (2.40) 

Here, the correlation length ^ ~ 1/AE is related to the lowest energy gap in the hamiltonian 
H, such that the dynamical exponent z = 1. 

Finally, for the correlation function, we have from fl2.27p that at criticality, where a = 0 


C(R) = {S„Sn) = (4^) B'-- (2.41) 

In contrast to the thermodynamics observables considered before, this result^ holds true for 
arbitrary dimensions and is not restricted to d < 3. 

^Observe that the exponents of R in C{R) ^ R~^‘^~^'> for ^ ^ i? and C{R) ~ for ^ <c i? are 

different. For d = 2, one recovers the Ornstein-Zernicke form. 


12 


























Table 1: Critical exponents for the quantum spherical model fll.lll) at zero temperature, along 
the quantum critical isochore i? = 0, in dependence on the dimension d and the coupling A. 


critical isochore 


a 

/S 

7 

V 

V 

7 

d < 3 A ^ 0 

(d- 

3)/(<i-l) 

1/2 

1 — 1 

1 

cT 

1 — 1 

1 

1 — 1 

0 

1 

d > 3 A ^ 0 


0 

1/2 

1 

1/2 

0 

1 

d > 0 A = 0 


0 

1/2 

1 

— 

— 

2 


Table 2: Critical exponents for the quantum spherical model (11.111) at zero temperature, along 
the quantum critical ‘isotherm’ tg = 0, in dependence on the dimension d and the coupling A. 


critical isotherm 


Otc 

7c 

d 

I'c 

d < 3 
d > 3 

o o 

2(d- 

3)/(d + 3) 

0 

4/(d + 3) 
2/3 

(d + 3)/(d-l) 

3 

CO —|— 

CO 

d > 0 

A = 0 


0 

-1/3 

3 

- 


For the interpretation of these results, we recall the conventional critical exponents and also 
the associated amplitudes, in the notation of [60] along the quantum critical ‘isochore’ B = 0 


A 


m 


^Dt^n ; X ^ rt/ ; c ~ -t “ + Co ; ^ ~ ^ot„ ; G{R) ~ ^ ; 


a 


The values of the exponents can be read off and are collected in table [TJ 
agree with those of the classical spherical model in d + 1 dimensions. 

Along the quantum critical isotherm, = 0, one can define 


AE ~ (2.42) 

As expected they 


c ~ Co + [Ac/ac)\B\ ; y ~ rc|i?| i? ~ Dcrn\m\^ ^ ~ ic\B\ (2.43) 

and read off the exponents]^ collected in table [2l The universality of this quantum phase 
transition is conhrmed through the A-independence of all these exponents. 

In addition, the universality of full scaling scaling forms fl2.26112.29p can be tested by working 
out at least three universal amplitude combinations [60] • Considering the singular free energy 
and its derivatives, we considered three amplitude combinations which from fl2.26p are expected 
to be universal. Explicitly 


= (1^ 

Rx = TDcB^-^ = 1 

6r,Dl/^ = 1 


(2.44) 


®In order to avoid ambiguities, we write D for the amplitude denoted as B in |60j . since we have already 
used the letter B to denote the magnetic field. Analogously, along the quantum critical ‘isotherm’ tg = 0, we 
write Dc instead of the conventional notation Be |60j . 

®These obey the standard scaling relations, such as Uc = a/(35, 7c = 1 — 1/^, r'c = v/(35. 
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and we give the results which follow from our explicit calculations above. The A-independence 
of these three amplitude ratios is additional conhrmation of the scaling form fl2.26l) . with only 
two non-universal metric factors. In order to test the universality of the scaling form fl2.29p of 
the spin-spin correlator, consider 


Qi 


Qs 






d+1 


r jf+(o) 


wyioy- 


(2.45) 


whose universality is conhrmed explicitly through the A-independence. Observe that for 1 < 
d < 3 all universal amplitude ratios in fl2.44112.45p are hnite, but that several of them they either 
vanish or explode when d —?■ 1 or d —)■ 3. This indicates that the scaling behaviour is going to 
be different (or does not even exist) when d > 3 or d < 1. 

For the spin-anisotropic quantum spherical model, we can conclude that the scaling forms 
! o.nd their universality, have been fully confirmed at the quantum critical point at 
T = 0, g = gc{X,d), with 1 < d < 3 and 0 < A < 1. Since the scaling functions themselves 
are universal, they were already calculated explicitly in the classical spherical model in d -|- 1 
dimensions, see e.g. m, and need not be repeated here. 

2. For 0 < A < 1 and d = 3, we are working at the upper critical dimension. Therefore, 
we have to introduce logarithmic corrections to the scaling behaviour, see eq. fl2.3ip . In order 
to work with the logarithmic terms and the magnetic held, we introduce the dimensionless 
held B := ^ I'his manner, the expression Ini? is well-dehned. We hnd for the 
magnetisation 


m{tg, 0 ) ~ 

m(0,i?) ~ 


8J 




i?| 31 In |i?| 13 


and for the susceptibility 

x(0,5) ~ ^ ■ \B\~^\n\B\\^ 

In the same manner as above, we calculate the specihc heat and hnd 

c(i3,0) ^ ■ \ In\tg\\-^ 

c(0,5) ~ 3^^^-I In 1511-^ 


(2.46) 

(2.47) 

(2.48) 

(2.49) 

(2.50) 

(2.51) 
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Finally, the correlation length reads 


■ l^sl (2.52) 

^(0,5) ~ \l~^ (^) ■ 

This logarithmic behaviour can be described in terms of logarithmic sub-scaling exponents [H] 

c ~ \tg\~°'\\n\tg\\°‘ ; m{tg,0) ^ \tgf\\n\tg\f ] X ^ \tg\~'*\'^^\tg\\^ ■, 

e ~ \tg\-'^\\n\tg\f ; m{0,B) B^/^\ln\B\f ; C{R) (2.54) 

and we simply read off their (universal, since A-independent) values 


« =-1 ; 0 ; 7 = 1 ; F= ^ ^ ; 77 = 0 (2.55) 

These values agree with those of the 4iA 0(?7,)-Heisenberg model in the limit n —?■ cxd [ISl ITT] . 

3. In the case 0 < A < 1 and 3 < d we expect mean-held critical behaviour. Near criticality 
0 < fg 1, we hnd the observables in the same manner as in the previous parts, but with the 
’linear’ spherical constraint. We hnd the observables along the critical 5 = 0 line 


m{tg,0) 

xitg,0) 

c{tg,0) 


atg,0) 




SJ 


■tg^ 


A 


2J 


^g 


V2A, 


-y^K^Jgc 



and along the quantum critical isotherm tg = 0 they read 


m(0, B) 

X(0,5) 

c(0,5) 

m.B) 


2JA- 


8J 




1 

' Al> 

1 8J ’ 

2J 

(2J)2\ 



2 

■ 5-3 


V2A, 




1 + A 


(2J)3h^ 


2Al. 


■ 5-3 


(2.56) 

(2.57) 

(2.58) 

(2.59) 

(2.60) 

(2.61) 

(2.62) 

(2.63) 


Reading oh the critical exponents (see tables [T] and [2]) yields the expected mean-held behaviour. 
4. For A = 0 and d arbitrary, the free energy density reads 


fitg,B) 


52 1 

TT-^ ^ 

4J a 



(2.64) 
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The magnetisation reads consequently 


m{0,B) ~ 

and the magnetic susceptibility becomes 

X{tg,0) ~ 

X{0,B) ~ 


i ■ t/'" (2,65) 

(2 ( 2 . 66 ) 


1 Vs , 

2J d ‘ 

(2.67) 


(2.68) 


The specific heat is found to be constant near criticality and along the quantum critical isotherm 


c ~ 



(2.69) 


The critical exponents are listed in tables [T] and [2J They are distinct from those of the modified 
quantum spherical models defined in [211132], where the particle number N is conserved as well. 

For the correlation function, we see a disconnected part from the zero temperature contri¬ 
bution. As derived in appendix D, we have to take thermal contributions into account. We 
then find 

^ ^ (-2^s) Io{zY-^Ir{z) (2.70) 

with z = gJh?/2T'^s. At criticality, we can deduce to leading order in T, see eq. (1D.17P 


C{R) 


AdJ 


<5r,0 + 


dJ [ 


, T^d 


T 


■2-d , 







(2.71) 


with the thermal reference length T'^d/gcJh^ and where the critical coupling constant 

Qc = gc{0,d]T) has to be found from the spherical constraint in the non-vanishing zero- 
temperature limit. To leading order in T, this gives the condition 



fT , / d 

y Ad (27r)'^/2 {g^JK^ J 


(2.72) 


hence gc — Ad{l — illustrates how finite-temperature effects 

renormalise the value of g^- The behaviour fl2.7ip of the correlation function does not fit into 
the standard phenomenology, described by the conventional critical exponents [211391 Ell H]. 
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2.5 Casimir effect in d = 1 dimension 


Although an analysis of finite-size effects is beyond the scope of this work, we add a brief 
comment on the Casimir effect in the d \ 1 limit, that is the case of a strip geometry, of finite 
width L, and with periodic boundary conditions. 

For ID quantum systems with sufficiently short-ranged interactions and a classical corre¬ 
spondent model such that z = 1, conformal invariance is expected to hold at the quantum 
ciritcal point at temperature T = 0, see [24l |39]. Scale-invariance alone gives for the nor¬ 
malised free energy density f/D q = fo — L~‘^Y+ o(L“^) where Y is 
an universal scaling function and Ci ^2 and Dq are the non-universal metric factors [59l SO] . The 
normalisation constant Dq must be fixed such that the dispersion (energy-momentum) relation 
becomes E{k) = \k\ for fc —)■ 0, snch that energy and momenta are measured in the same units, 
see [39]. Then conformal invariance relates the universal value y(0,0) = —7rc/6 to the central 
charge of the corresponding 2D conformal field-theory |2]. For the quantum XY chain (II.Sh . 
f /Dq has indeed been calculated, Y(0,0) was shown to be universal and the central charge 
c = I was found [37|, as expected for a model in the universality class of the 2D classical Ising 
model [21I391125]. 

If we want to apply the same method to the quantum spherical model in d \ 1 dimensions, 
we have to take into account the possibility that the critical value Sc of the spherical param¬ 
eter may acquire a finite-size correction. Explicit calculations have shown, however, that this 
universal finite-size amplitnde vanishes, for periodic bonndary conditions, when d \ 1 [siiin]. 
Hence, f /Dq can be taken over from the free fermion representation of the quantnm XY chain, 
where the boson-fermion correspondence implies that periodic boundary condition in the even 
sector (to which the ground-state belongs) of the quantum spherical model corresponds to 
anfFperiodic boundary conditions in the even sector of the fermionic model fll.Qp |53] . Hence, 
the ground-state energy of the periodic spherical model chain is identical to the ground-state 
energy of the quantum XY chain, with anti-periodic boundary conditions. This is known to 
read [371 EH] 

= + + (2,73) 

where /o(A) is an explicitly known, non-universal bulk contribution to the free energy density. 
We see that the finite-size amplitude Y(0,0) = —vr/b is A-independent and therefore universal, 
as expected [59] HO] , but the higher-order finite-size corrections are non-universal. We find the 
value c = 1 for the central charge in d \ 1 dimensions, as expected for a free boson. 

For dimensions d > 1, the simplifications we could use here, in the d \ 1 limit, do no 
longer apply such that the computation of the Casimir effect is considerably more involved, see 
[MumiHiiiiiiiiniiisiiHoiiii] and references therein. It would be interesting if recent attempts 
to formulate a conformal bootstrap for the 3D Ising model [26] could be brought to shed light 
on the interpretation of universal Casimir amplitudes. 


3 Conclusions 

We have explored the T = 0 quantum critical behaviour of the spin-anisotropic qnantnm 
spherical model fll.lip . One of our motivations was to be able to compare the effects of bosonic 
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versus fermionic degrees of freedom, by using the information available from the quantum 
XY model [IHl El EH IMl EH EH EH EH ES]. However, the quantum spherical model has the 
advantage that it can be analysed exactly for arbitrary dimensions d, coupling g and external 
held B, whereas the quantum XY model is only solved for d = 1 and for a vanishing external 
held B = 0. As to be expected, we have found a line (‘quantum critical isochore’) of quantum 
phase transitions and used the pair creation/annihilation rate A > 0 to test explicitly for 
universality along this line. The generalised Privman-Fisher scaling form, adapted to quantum 
criticality [59l EOl El] allowed to test not only the universality of the exponents but also of 
certain universal amplitude rations and in consequence of the full scaling forms fl2.25ll2.27p . It 
is known since a long time that the critical behaviour of the fermionic model along the critical 
isochore is universal [361 EH; we obtained here the analogous result for the bosonic model. 
Merely the values of the exponents are diherent (in ID, the identihed central charges also 
diher). In the quantum spherical model, an analogous test can also be carried out along the 
quantum critical isotherm tg = 0. 

In the special case A = 0, the total particle number is conserved, leading to a different global 
symmetry and the critical behaviour is different. It is also distinct from the spherical model 
variants [56l 02] with a global conservation of the number of quantum particles. 

In the fermionic quantum XY model, the ordered phase contains a sub-phase, for 0 < ^ < 
d\/l — A^ , with spatially oscillating correlation functions [6] EH EH] EH EH- This sub-phase 
is characterised by level crossings in the hamiltonian energy spectrum, between the even and 
odd spin sectors |12], The transition line between oscillating and non-oscillating correlators, at 
g = dVl - A2 , is characterised by the existence of certain Neel ground-states [5T]- We did not 
succeed to detect similar properties in the bosonic quantum spherical model. 

A surprising feature of the model studied here is the re-entrant quantum phase transition 
for dimensions d < 2.065 and sufficiently small values of A. This shape of the quantum critical 
line could not have been anticipated from previous studies of the classical spherical model. This 
makes it clear that interactions between the momenta cannot always be absorbed into a change 
of variables^ In hgureE] we compare the shape of the critical line gc = 5'c(A), normalised to 
the value at 5'c(0) at A = 0, of the bosonic quantum spherical model fll.lip . with the fermionic 
quantum XY model. In ID, the latter model reduces to free fermions. Comparing the shapes 
of gc{X,d), the re-entrant phase transition found in the bosonic case of the SAQSM does not 
appear in the analogous ID fermionic model, where gd^) = 1 is simply constant [H EH- 
order to better appreciate the influence of dimensionality in the quantum XY chain on fi'c(A), 
and in the absence of an analytic solution, the best what we can do is to compare with the few 
known numerical values of gd^) in extension of the spin hamiltonian Hxy from fll.Sp to 2D 
[36] . Although those few data shown in £gure[3]seem to indicate that the approach of gd^) 
towards the A = 0 case should be monotonous and hence no re-entrant transition is suggested, 
the available data are too few and too far apart for a final conclusion. 

Since in many respects, effectively non-integer values of the dimension d can also be pro¬ 
duced by long-ranged interactions [SSI EH EH EH 12B] , one could anticipate that several of our 
conclusions might have qualitative analogues in long-ranged quantum phase transitions. Also, 
it would be interesting to see if the theory of random matrices, so sucessfully used in fermionic 

Considering the leading finite-temperature corrections to the value of gc(A, d), it can be shown that for T 
sufficiently small, the value of Qc is only slightly renormalised such that the re-entrant transition also occurs for 
finite (and small) temperatures T > 0. 
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Figure 3: Left panel: normalised critical coupling gd,^)/QciS^) in the quantum XY model fll.Sp . 
as a function of the coupling A. In ID, one has = 1. In 2D, the numerically known esti¬ 
mates of gd^) [3n] are given by the dots and the dashed line is a guide to the eye. Right panel: 
normalised critical coupling gd^,d)/gdO,d) in the quantum spherical model fll.lip . as a func¬ 
tion of A and for dimensions d = [1.3,1.4,1.5,1.6,1.7,1.8,1.9, 2.0, 2.1, 2.5, 3.0] from bottom to 
top. 


quantum chains [31111], could be brought to be applied to the kind of bosonic systems analysed 
here. 

This illustrates that the interactions between the conjugate momenta can play a physically 
important role. Our results raise the question of the quantitative importance of more general 
kinetic terms, e.g. in 0(u)-symmetric quantum rotor models with n hnite. Also, one may 
anticipate a rich phenomenology when combining different kinds of interactions between the 
spins and the momenta. If such effects should be found, the spherical model would have 
demonstrated once more its usefulness as a heuristic device and guide towards non-trivial and 
interesting new types of critical behaviour. 


19 




























Appendix A. Diagonalisation via a canonical transforma¬ 
tion 


The quantum hamiltonians to be diagonalised are of the form 






+ ^hBnmO^ln) 


^ ^ C'n (On + ®|i) 


(A.l) 


where the sums run over the J\f = N'^ sites of a d-dimensional hyper-cubic lattice. A is 
a hermitian matrix, B a symmetric matrix and C is a real, constant vector. The bosonic 
annihilation and creation operators Un, obey the standard commutators 

[cZfj, Qim] ^ ^nm (A.2) 

For (7 = 0, the diagonalisation procedure follows closely the fermionic techniques of Lieb, 
Schultz and Mattis [53], applied to quantum Ising/XY chains. In the bosonic case, any space 
dimension d can be treated and (7 7 ^ 0 is admissible. Throughout, we restrict to the case when 
A and B are real-valued (although extensions are readily formulated). 

We seek a canonical transformation which brings H to the form 

H = Afc + Hq (A.3) 

where are again bosonic annihilation/creation operators, are the sought eigenvalues 

and the constant Hq has to be determined. The required canonical transformation is of the 
form 

bk '^k T ^ ^ (VkpCip Wkptt^pj , bf^ Uk T ^ ^ T (A. 4 ) 

p p 

where the Af x J\f matrices V and W are determined from the bosonic commutation relations 
and the Uk are numbers. This gives VV^ — WW^ = and VW'^ — WV'^ = 0, where ^ 
denotes the transpose and is the N'^ x unit matrix. A direct consequence of these is 
iy + W){V- Wf = Id, hence 

{y + w)-^ = {y-wY , {y-w)-^ = {y+ wf (A. 5 ) 

The last conditions on V, W come from the requirement that the canonical transformation flA.4p 
brings II to its diagonal form flA.3p . which means [5fc,i7] = A^bk- Hence 

y±w)(A±B) = AyTiv) , y-w)c = Au (A.e) 

where A = diag(Ai,..., A^y) is a diagonal matrix with the eigenvalues A^, and the vector 
u = (ui...,u^). Following [53], one dehnes two matrices, arranged as two sets of vectors 
($fc)m := y + W)kni and := {V - hF)fcm so that by reading eq. (IA. 6 D line by line, one 

has the two coupled equations 

iA + B)= Ak<I^l , <IfliA-B)= Ak^l (A. 7 ) 
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so that the eigenvalues can be found from the following eigenvalue equatioiJ^ 


:= (A -B){A + B) = AlA^l 


(A.8) 


Later on, we shall also need the explicit transformation of the creation/annihilation operators. 
For Cn = 0, this reads b = Va + WTA and its inverse becomes a = V'^b — W'^b\ along with the 
hermitian conjugates. We shall require this below for the calculation of correlators. 


Next, we must hnd the eigenvalues A^ for the specific hamiltonian fll.lip in the main text, 
with nearest-neighbour interaction^^. Then the diagonal form of the hamiltonian is 

given by hA.A\) . where the eigenvalues A^ = sKk are, for a hyper-cubic sguare of M = N'^ sites 
in d spatial dimensions and with periodic boundary conditions 


Afc — 



1 + A 
2 s 


d 

cos kj 
i=i 



1-A 
2 s 


cos kj 




(A.9) 


where the guasi-momenta kj = -§-nj, with Uj = 0,1,... N — 1 and j = 1,... ,d and the spherical 
parameter s. 


Proof: Eq. flA.9p can be derived from the properties of cyclic matrices |1] and using mathe¬ 
matical induction over the dimension d. In what follows, we denote a cyclic Af x J\f matrix, 
generated from a vector (t>i,..., by 


^{vi, ...,vm) 


Vi 

V2 

Vs ■' 

■ ■ VM-1 

VM 

VM 

Vl 

V2 ■ ■ 

■ ■ Vj,f_2 

VM-l 

V2 

V3 

V4 ■ ■ 

■ ■ Vj^ 

Vl 


For the sake of this proof, we work with the reduced, dimensionless hamiltonian Hr := 

H/ih^fgjT). 

Step 1 : For d = 1, M = N. The matrices A = A^^'> and B = are (the index refers to 
the value of d) 


4s’ 


4s 


A« = e: i,--,o,...,o,-- , = o,-,o,...,o,- 


A 


A 


4s' 


4s 


and therefore 




8 s2 ’ 2s’ 16s2 


16s2 ’ 2s 


is cyclic as well [Ij. The eigenvalue equation flA.8p can now be solved by the ansatz (4/^)^ = 
Since the cyclicity of all matrices implies periodic boundary conditions, this produces flA.9D for 

^^The only difference with respect to fermionic chains [53] is that therein B = —B^ is antisymmetric. For 
bosonic as well as for fermionic systems, the matrix {A — B)[A-\- B) is symmetric and positive semi-definite, 
such that all eigenvalues are real. 

^^The method outlined in this appendix works for arbitrary interactions, although the practical calculations 
can become more involved. 
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d = 1 and the values of k are indicated!^ 

Step 2: In order to demonstrate the passage from dtod+1 dimensions, consider a multi-index 
notation in d -|- 1 dimensions 

n={ni,n 2 ,...,nd,nd+i) = innd+i) , n = (ni, ua,..., n^) 

where individually, nj = 0,1,..., — 1, with j = 1, 2,..., d, d -|- 1. In d -|- 1 dimensions, the 

hamiltonian can be brought to a block form as follows 




d+l 


J = 1 


^ ^ I yOi-nO^n+ej + ^n^n+Sj j “1“ ®li®n+ej + dnCLn+Sj 


N-1 


EE 


k=0 

d+l 


(^lik^nk + 


E/ {^nkO-nk+ejk + + o!~ik^nk+ejk + O^nkO^nk+ejk 

t=l 


4s - 

Af -1 


{<^nk<^nk+l + 


E/ "i ^ Ao E/ I + <^\ik^^nk+^ + + ^■nkCi'^k+l 


k=0 


N-1 


n,m k^i=0 


As 


)n^ _ - n~ _ 

/ . / . ] ^nk^nk,me^rne ^ ^rikJ^nk^me^'mi “T ^nk^Hk,me^7fU 


(d+l) 




2 L 


where Hjf^ is the local hamiltonian in the fc-th d-dimensional layer. The interaction matrices 


have the block structure 


4tS, = 444/- 

4114 = 444/ + ^ (4/+1 + 4/-i) 1/ 

where is the iV'^ x unit matrix. In turn, they may be written as cyclic matrices of blocks 


4s’ 


4s 


= ^ ^(4,_ ,o,...,o,-— , = ,o,...,o,— 


A 


A 


2s’ 


2s 


Next, we write down the block structure of the eigenvalue equation (1A.8I1 


.M+i) + 1-A2 l-A^ + XB^'^^X 

^(rf+i) ^ ^(4 +---^ 0 ,..., 0 ,-- 

\ 8s 2s Ibs^ 16s^ 2s / 

Now, the habitual ansatz = ^fi where by induction hypothesis, is the eigenvector 
of the d-dimensional problem, gives for the eigenvalue in d -|- 1 dimensions 




= + 


2 1 -A^ 


4s2 


cos^ kd+i 


cos k, 


d+l 


1 _ A 2 “ 

1 -::-cos kj 


2s 


i=i 


13 


For kj with j 4 0, iV/2, the eigenvalues = Xkpf_j are degenerate, so that the corresponding eigenvectors 
can always be chosen with real-valued components. 
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COS k, 


d+l 


1 + A 

1-;— y. cos kj 


2s 


1 — A / 1 — A \ 

+ I 1-— 2. cos kj 


1 + A 


1 1 + ^^ . 

1 -;— y .cos k 


2s 


i=i 



2 

d+l 


2 s 


i=i 


1-cos /Cj 


2 s 


i=i 


where in the last step the induction hypothesis flA.Op was used for in d dimensions. This 


completes the proof. q.e.d. 

Finally, we find the constant Hq in (IA.3F For the sake of notational simplicity, we only treat 
the case d = 1 explicitly, but we shall give the generic result at the end. Since the eigenvalues 
are generically two-fold degenerate, we first go over to real-valued combinations 

•= / l<^k[{'^k)n +i^Dn] = Ck cosnk ; if fc < A^/2 

\ jiCk[{'^k)n-i'^*k)n] = Ck sin nk ; if fc > A^/2 

Here, Ck is a constant which will provide appropriate normalisation. 

From (lA. 6 p . we further have 4*^ = A^^(A — B)'^k, hence 


(A.IO) 


(^^)n CfcA^ I 1 


1 + A 
2 s 


cos k 


cosnk ; if k < N/2 
sin n/c ; if A; > N/2 


(A.ll) 


The normalisation constants follow from the bosonic commutator relations and which require 

E + -+" = E (**)»(«'*)» = ! (A.12) 


so that finally 


At 


Cu = 


1/N ; if fc = 0, N/2 

1 — (1 + A)(2s)“^ cos A; \ 2/N ; else 


(A.13) 


The extension to d > 1 dimensions is now obvious. 


While this gives the general method, we now apply it to the specific hamiltonian fll.lip in 
the main text. For a spatially constant magnetic field, all constants are equal Cn = C. From 
eq. (IA. 6 p . we deduce 


C 


N-l 


^ n =0 


kin 


(A.14) 


Using the geometric sum, it is obvious that Uk vanishes for A; 7 ^ 0. For A; = 0, we find 


uq = CNKq^cq (A. 15) 

Thus, we are now able to write down the constant Hq by rewriting the diagonal hamiltonian 
in the form j + Ho-, using the transformation formula flA.4p and 

comparing the constant terms. We find Hq = —AqMo hence the ground-state energy reads 


Un = 


-aqmo+o 


C^N 


k — 


(1 + A)(2s) 


-1 




(A.16) 


with k = and n = 0,1,... iV — 1. The generalisation of flA.lbp to d > 1 is obvious. 
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Appendix B. Spherical constraint for 


We derive the spherical constraint eqs. fl2.12112.13p in the main text, for general A 7 ^ 0,1. 

Since the magnetic held term in ( 12 .9p is just additive, we can set i? = 0 for our purpose. 

Starting from the form (12. 9 p of the spherical constraint, the product of the two square roots 
in the denominator is folded into a single factor by the Feynman identity, see e.g. [5] 




1 

TT 


da: 


0 ■^x(l — x) xA + (1 — x)B 


(B.l) 


so that the constraint becomes (with the Brillouin zone B = [—tt, vr]'^ 

da: 


'Svr^J _3 

= s 2 


dfc 


2 _ 1-A2 

A 


cos kj 


h?g 


x{s-^ ZU + (1 - ^) (s - ^ EU cos kj') 

(B.2) 

However, we are looking for a representation which factorises in the momenta kj, such that the 
dimension d can be treated as a real parameter, in analogy to the known representations valid for 
A = 1. The denominator could be simply exponentiated, via the identity G~^ = J^due~^^, 
but the terms in the numerator still couple the different kj. One might consider to obtain 
these factors by deriving the exponential with respect to a: or 1 — a:, but this cannot be done 
immediately, since the presence of x in both terms in the exponential would generate unwanted 
contributions. It is better to introduce hrst an auxiliary variable 

y = 1 — X (B.3) 

and to render it formally independent from x, by inserting a Delta function into an additional 
integration over y, according to 

f dy 6{y — 1 + x)f{x,y) = f{x,l — x) , for 0 < a: < 1. (B.4) 

Jo 

Now, changing the order of integrations, we can indeed re-write the denominator as an expo¬ 
nential and afterwards express the numerator as a derivative of this exponential. This is done 
by dehning the differential operator 

A 


' udx udy 


su- 


dxdy ) 


(B.5) 


Then eq. flB.2p can be re-written as follows 
'SttV 


n?g 



\ 5{y + x-l) 

du — -^ X 




cos ki 


uy 


i=i 



xV. 


xy 


VW 

r / \ 1 T f ITA 1 — A 

exp [-u{x -F y)s\ Iq ux— -h uy—— 


(B. 6 ) 
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since now the integrations of the kj factorise and can be carried out separately. Here and below, 
the /n(p) denote modified Bessel functions [1]. 


Here, a further comment is necessary concerning the argument of the modified Bessel func¬ 
tion. Clearly, and taking into account that y = 1 — x will have to be put back, the argument 
vanishes linearly at 

xo = ^{l- A"^) . (B.7) 

For 0 < A < 1, one has Xq < 0 which is outside the interval of integration and need not concern 
us. But for A > 1, one would have 0 < xq < | inside the integration interval of x, since the 
derivatives of Jq lead to higher order modified Bessel functions with n > 1, which vanish for 
a vanishing argument. Then a more careful distinction of cases which takes these zeroes into 
account will become necessary. 


We now apply the operator V^y 


to the integrand in flB.6j) and also define 


g := g{u, x, X) = u 



X 


A 


(B.8) 


Then the spherical constraint flB.6j) becomes 




exp [—ms] 
\/x{l - x) 


[h{Q)Y X 


X 


1 — A^ 



dl-\^ 

2~^ 



(B.9) 


Eqs. flB.9|) and flB.8|) are eqs. (12.12112.131) in the main text. 


Indeed, the dimension d can now be considered as a real parameter, which offers obvious 
conceptual advantages. For s > Sc = ^^d and A 7^ 0, this integral is convergent for all d > 1. 
While this representation, as it stands, holds true for all values of A, the asymptotic analysis 
will become more simple for 0 < A < 1, where the possibility of zeroes of the In{g), with ri > 1, 
need not be taken into account. 


Appendix C. Asymptotic behaviour 


We analyse the spherical constraint fl2.12p and derive the asymptotic relations fl2.18p for generic 
couplings 0 < A < 1. 

1 . For 1 < d < 3, the leading contribution to the shift tg in the coupling g is non-analytic. 
Considering the spherical constraint (12.121) . non-analytic contributions come from large values 
of u in one of the integrals. Combining eqs. (12.12112.16p . we must analyse 


t 


9 



V\/d 


1 

^g3/2 



exp [—ms] 
\/x{l - x) 


[Hq)Y X 


X 


1 - A^ 

d(d-l)^ 



^ 1-A^ A , h{g)\ 
2 4 ^ ^ Ug)) 


(C.l) 
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for a spherical parameter s = Sc + <^=|(l + ^)d + a in the vicinity of a —?■ 0. Here 77 is a 
cut-off which helps to isolate the non-analytic contributions to tg and we shall let 77 —)■ cxo at 
the end. Because of fl2.13p . the argument q of the modified Bessel functions never vanishes for 
0 < A < 1. Then, in order to obtain the leading behaviour in a, it is enough to use the leading 
asymptotic behaviour /^(p) — /^2tiq (1 -|- 0(l/p)) [1] of the modihed Bessel functions. Then 

Inio)/I o{q) — 1 to leading order in 1 /p for n = 1,2 and we arrive at 


f 

Lg — 


1 


TT 


du / dx 


exp [—ms] f exp g A ' 


0 y'x(l - x) \V27rgJ 


s^ — d{d — 1 ) 


l-A' 


d 


l-A' 


For convenience,we recall the definition of g from fl2.13p 

1 +A 


g = p(m, X, X) = u 


X- 


+ (1 - x)- 


A 


2 ^ '2 

and absorb into a single constant k several purely numerical factors 

A 2 


K := 


d^ 


4 


—- — 1 


TT 2 


,3/2 


dA 2 


2(1 A) 7 r ^+'^/2 


(C.2) 


(C.3) 


(C.4) 


such that the constraint becomes more compactly 

, exp [—ms] exp lux^d] exp [m(1 — x)^d] 
da; 


ig — 


K 


du 


'V 


'0 


»1 


= K 


du / da; 


\/x{l — x) (1 — A(1 — 2xyf^’^ 

exp [—ua] exp [—m(1 — x)dX] 


0 ^a;(l - x) m '^/2 (1 - A(1 - 2x)) 


d/2 


= 2k duu ‘^1’^ exp [—ua] / d(j) 


2Kai-^ 


r-7r/2 


dip 


<0 (1 -|- A cos2ip) 


d/2 


exp [—udX cos^ (f 
(1 — A cos 20) 2 

CX) 

dn exp 


' rja 


vdX 


sin^ 0 


a 


(C.5) 


where we used s = |(1 -|- A)(a; + 1 — x)d -|- a in the 2 "^^ line and changed variables several times. 


in the 3™ line according to a; = sin^ 0 , and in the 4*“ line v = ua and 0 = f — 0 , and also used 
cos 0 = sin 0 and cos 20 = — cos 20 . 

We are interested in the asymptotic behaviour near criticality, when 0 < a 1. Further¬ 
more the main contribution to the 0 -integral, for v still hnite, will come from the region where 
sin^0/(T = 0(1). But in the a —?■ 0^ limit we consider here 0 will be small as well so that we 
can replace sin0 ~ 0. Then the main contribution to this particular integral in (IC.SP should 
come from the region 

02 < a . (C. 6 ) 

Hence the leading term can be obtained by replacing the upper limit in the 0-integral in flC.Sp 
by infinity. Changing the order of integrations, we find 


t„ ~ 2Ka2 


-1 


' rja 


dn 

yd/2 


exp 


d0 


(1 + A) 


d/2 


exp 


-vdX^ 


a 



^(d-i)/2 r (V,r7a) ^/X/2 
[77(1 -|- A)]C+i )/2 


(C.7) 
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with the incomplete Gamma function r(a,a:) pQ. Next, we have to carry out the two limiting 
processes, first a —)■ 0’*' and then r; —)■ cx), in exactly this order. Dehning the Gamma function 
via analytical continuation, for 1 < d < 3 we simply have linio-^o r(^^, a) = r(i^) and obtain 


, ^ ,(.-!)/. r (1^) v/a72 

[7r(l + A)]G+b/2 


Al< for 1 < d < 3 


(G. 8 ) 


2. For d = 3, we can repeat the analysis leading to flG.7p . However, the limit a —?■ O’*" in the 
incomplete Gamma function has to be taken more carefully. Using [U eqs. (6.5.19, 5.1.11)], 
one has a logarithmic term 

r(-l, x) ~ - + Gb - 1 + ln(a:) - ^ + 0(a;^) (G.9) 

X 2 

where Ce ~ 0.5772... is Euler’s constant. Gonsequently, we find for the a-dependence in tg 

(jr(— 1 , rja) ~ —h [Ce — 1 + ln(? 7 (T)] a + O(cr^) ~ cxlncr (G.IO) 

In the last expression, we merely retain the most singular term when a —)■ 0^ with rj finite 
and then dropped those terms which vanish in the r] ^ oo limit. The leading non-analytic 
contribution in flG.7ll is 


f 

Lg — 


yv2 


|ir(l + A)]; 


(T In (T =: H 3 (T In (T ; for d = 3 


(G.ll) 


3. For d > 3, the non-analytic contribution from eq. flG.7p tg ~ is of higher 

order than linear. The leading term in tg now comes from the the analytic contributions to 
fl2.12p which was previously subtracted from the left-hand side. The leading correction term 
is found by a straightforward expansion in a. We also introduce the short-hand F{d,X,p) := 

—d(d — ~ obviously independent of s. Hence, recalling 

also flG.3p 



[ du [ o(g) _ ^ 3/2 [—ms] (s^ -|- F(d, A, g)) 

0 Jo ^/x{l-x) 

0 do \/x{l-x) 










In this expansion, the zeroth order gives Qc and the first order gives the required linear contri¬ 
bution tg ~ Ay a, where H> is given below in flG.15p . Its value must be found numerically. 

Summarising, we have found, for 0 < A < 1 


t 


9 




H 3 cr In cr 
Ay a 


if 1 < d < 3 
if d = 3 
if d > 3 


(G.12) 
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with the following constant amplitudes (derived here for 0 < A < 1 but which can be continued 
to A = 1 as well) 


Al< 


^3 




r(¥) 

[ 7 r(l + A )]('^+^)/2 

A 72 

l7r(l + AP 

^ Jo Jo V^(l - x) 


2uSc — 1 


3 2xiSq 


F{d, A, g) 


e-us. 


(C.13) 

(C.14) 

(C.15) 


with Sc = {I + X)d/2, F{d, X, p) was defined above and flC.3p was used. On the other hand, for 
A > 1 the argument g of the In{g), as given by (1C.3I1 . can vanish, the analysis leading to (1C.7I1 
has to be re-done and (10.12^ cannot be expected to remain valid. 


Appendix D. Spin-spin correlator 


Using the representation fll.Sp in terms of ladder operators and then the canonical transforma¬ 
tion flA.4IIA.8p from appendix A, the spin-spin correlator is given by 


{Sr^S^) 



(D.l) 

(D.2) 


Since the ladder operators are bosonic, they obey Bose-Einstein-statistics. Hence 
(bkh'^ = = 0 ; = 6k,k' (exp [Afc/T] - 1)”^ 

This immediately leads to 


{Sr^Sm) 



^(^'fc)r.(^fc)m COth [Afc/(2T)] 


(D.3) 


(D.4) 


Using the real representation of the vector from appendix A, we find for the correlator in 
the continuum limit, with m = n + r 


{SnSn+r) — 


h?g 


dk 


2s - (1 - A) cos kj 


—/ 7 — 7:7 A - 7 -coth [Afc/( 2 T)] TT cos (r,fc,) (D.5) 

SJsJe(27rr\2s-(l + X)j:‘‘,,coskj ' ” ''hi ' > 


and spatial translation-invariance is explicit, so that we can set n = 0 from now on. Eq. flD.Sp 
is an exact expression for any temperature T. 

1 . For A 7 ^ 0, consider the quantum phase transition at T = 0. Then fID.Sp simplihes to 


(SoSr) = 




dk 


2s - {I - X)Yfj^^cosk^ A 

[[cos [rjkj) 


8J s Ji 3 {2ttY ^ 2s - (1 -1- A) Yfj=i cos kj 


(D.6) 
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Because of explicit rotation-invariance, we can choose axes such that r = {R, 0 ..., 0). Now, 
eq. (ID. 6 P can be factorised by the same techniques as used in appendix B to factorise the 
spherical constraint. We hnd, with g from eq. flC.3j) 


(SoSr) = 


h?g 1 
8Js TT 


du 


dx exp [—ns] 

\/a;(l - x) _ 


1-A 


(d-1) 


h{Q) I'r{q) 

Io{q) Ir{q)) _ 




(D.7) 

In order to work out the correlator from this representation, we now analyse the main 
contributions to the n-integral. Since the integrand vanishes for n = 0 and u = oo, it will have 
a maximum at some intermediate value Umax and if the integrand is sufficiently peaked around 
Wmax; this will give the main contribution. Now, the leading term of the series expansion 

shows that for u not too large, the 
R. Since we merely interested in the 


Ir{p) ~ (p/2)^/r(i? -I- 1), for small arguments p <C 1, 
integrand will roughly behave as u^e~'^ such that Umax ~ 


large-i? limit, it follows that the contribution of small values of u to the integral is negligible to 
leading order. Therefore, in order to estimate {SqSr), it is enough to use the p 1 asymptotic 


form Rip) — (Svrp) ^/^exp 


P- 


of the Bessel functions, such that 






du / dx . = 

0 Jo V^(l-x) 


exp 


ag 


This can be evaluated following the lines of appendix C. We hnd 

-(d+l)/2 


IQ Q \ I ^^3 (1 + ^) (s (d-l )/2 I ^ 

This equation can be rewritten, using the naenuiuj- ) — 2 \ 2 ) 
for the modihed Bessel function of the second kind, to obtai 

^ 2-^!^ s-^d /^\('^-b/2 


-V — 


2g 


R^a 1 




(1 + A) 


V 


(D.9) 


the identity KRx) = A (f dvv ^ ^ exp [—v — x"^ j (4n)] 
he second kind, to obtain^q 




Ml 


(D.IO) 


and where the correlation length was identihed as ^ := Very close to criticality, ^ 

diverges, hence RR ^1. At some hnite distance from pc, one has on the contrary RR ^ 1 . 
Now, using the leading expansions [H eqs. (9. 6 .9,9.7.2)], one has the asymptotic behaviour 


TW J 2-^!^ T {i^) ■ R'-‘ ; 

' “ ™ V Js + \)d I 2 -'*''^^^e-*« ; 

(D.ll) 

with s = Sc + c, a is related to ^ and the value of g has to be taken from the spherical constraint. 

2 . If A = 0, we have to do a more careful analysis, since the zero-temperature contribution 
is completely disconnected. From eq. fID.bp . we see a ^r^q contribution arising. Thus, the 
leading non-trivial contributions in this particular case are thermal and we have to re-investigate 
the correlation function for non-zero temperatures. Hence we return to eq. fID.Sp . as well as 

^^For A = 1, eq. (ID.101) reproduces the well-known result [58l eq.(13)], if one takes into account that because of 
the normalisation RRn Ri) ~ W/4 chosen in |58] . one must renormalise s 1 —>• s/4 to ensure matching pre-factors. 
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to fl2.5p for the spherical constraint, in order to hnd the thermal corrections to the critical 
coupling constant Qc- Since we are still interested in a certain low-temperature limit and not 
in the thermal transition, we take 0 < T 1 and use the asymptotic expansion cothx ~ 
1 -|- 2exp(—2a:) to obtain the leading correction. The spherical constraint in zero held then 
reads _ 


1 = 


hg 


+ 2 


h'^g 


exp {—2zs) Jo {zY 


1 + 


d Ii {z) 
2s Jo {z) 


(D.12) 


8Js V 8Js 

with the argument := gJh?/2T‘^s. In the low-temperature limit, z ^ oo. From the 

asymptotic expansion of the modihed Bessel functions, we hnd 



2 e- 2 <"^ 

(2'kzY/‘^ 




(D.13) 


Studying this equation up to the leading order in 1/z, at the quantum critical point a = 0, we 
deduce the implicit equation for the critical coupling constant g^ = gd^, d; T) 



[T f d 

V 4(i (2ir)|'/2 ygjtfi j 


(D.14) 


First of all we see, that this equation is consistent with the zero-temperature limit and repro¬ 
duces gdO, d; 0) = AdJjtd correctly. While for d = 2, there is a simple closed solution 

9c{0,2-T) = 

eq. (1D.14P cannot be solved in closed form in general. 






(D.15) 


For large distances, the same techniques as before, applied to fID.Sp . lead for A = 0 to 


{SoSr) = \I^ + 


h?g 


exp {-2zs) Io{zY ^Ir\ 


z) 


(D.16) 


8Js V 2Js 

Using the asymptotic expansion for the Bessel functions, we hnd at the critical point g = g, 

d/4 


/c c \ _ , , jh^gc ( T^d Y 

“ V 457 + V17 iilTTfiS ) 


exp 


R^T 


d 


gcJf9 


(D.17) 


Appendix E. Critical coupling d) close to A = 0 


In order to prove fl2.24p and to understand the unexpected behaviour of the function gd^-, d) 
close to A = 0, we re-investigate the equation (recall the dehnition fIB.Sp of p = p(m. A, d)) 


J7r^(l -|- XYd'^ 

t9gd\ d) 


t t+A 


'o Jo Yx{l-x) 




1 +A 


-d{d — 1) 


1 - A2 f ldg) 
do{g) 


dl-X^ f IdQ) 
2 4 ^ ^ IdQ) 


(E.l) 


J'lr'^d^ 

hIgdO, d) 


+ 
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where the two contributions and G*^^) describe the leading behaviour in A, which are non- 
analytic and analytic, respectively. 


First, we consider the case 1 < d < 2, when the leading behaviour is given by the non- 
analytic term G*'^^ After a change of variable g = | (1 — A + 2Aa;) in flE.ip . we divide the 
^-integral in two parts d^ = d^) + d^?. In the limit A —)■ O'*’ and rj ^ oo, the first 

integral reduces to gdO, d) while the second integral will give the desired non-analytic term 
G*'^^ for small A. As in appendix C, G^^^ is analysed via the asymptotic expansions of the 
Bessel functions [1], which gives 


g(i) = 




dx exp [-2dgXj^^] 
'o \/x{l - x) (1 - A -h 2Aa;)(27r^)'^/2 


d^\ 


= 2 d 2 A '^/2 


t/2 


exp [—2dy sird'd~\ 
(27r|/)'^/2 


= 2d2A"/M dyexp[-dy]Io{yd){2ny)-^/^ 


(I'd / dy 
J Xt] 


(E.2) 


where in the second line, we made the substitutions y = Xg and x = cost'd and in the third 
line recalled the identity sin^ d = |(1 — cos 2'd) to derive dd = | e“"^/^/o(A/2) from 

the defining integral representation of Io{x) pp. For 1 < d < 2, this is indeed the leading 
contribution. Explicitly, using [611 cq. (2.15.3.3)], this further simplihes to 


fi 


gc{X,d) 


1/2^ J_. r(i-rf/2)r((rf-i)/2) 

Vid 2ir(‘i+i)/2 r(o!/2) 


(E.3) 


with a finite, positive amplitude for all dimensions 1 < d < 2. 

For d > 2, the non-analytic contribution G^^^ in flE.3F analytically continued in d, is dom¬ 
inated by a new analytic contribution To obtain this, one must formally expand the 

integrand in flE.lIl to first order in A. Of course, such as a formal expansion is only admissible 
up to the order where the expansion coefficient(s) converge(s). Because of the definition fIB.Sp 
of g, in principle the Bessel functions In{g) should be expanded around A = 0. However, the 
leading term will introduce a factor 1 — 2a; into the integrand and all these contributions vanish 
because of da; (1 — 2a;)/y^a:(l — x) = 0. Therefore, the additional contribution reads 


(i - #)) + + o(V) 

\ ^0(2)/ 2 ^0(2) J 

(E.4) 

and the integral over x has become trivial in the A —)■ 0 limit. This contribution is linear in 
A and hence will dominate over G^^^ for d > 2. In order to study its convergence, we split as 
usual ff^dg = f^dg + f°^dg and analyse the convergence of the second integral. Using the 
asymptotic expansion of the In(g) up to next-to-leading term in 1/g [1], the large-r^ behaviour 
of G^^^ is given by — ^(27r)“'^/^and this converges for d > 2. For d < 2 however, 
the integral diverges such that the formal expansion used to derive it does not exist. Then 
(IE.3D gives indeed the leading contribution to gc{X,d) for A 1. 

This proves fl2.24p in the main text. 
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